On the Expected Value of the Minimum 

Assignment 



Marshall W. Buck Clara S. Chan David P. Robbins * 

May 10, 2000 



Abstract 

The minimum /c-assignment of an m x n matrix X is the minimum 
sum of k entries of X, no two of which belong to the same row or col- 
umn. If X is generated by choosing each entry independently from the 
exponential distribution with mean 1, then Coppersmith and Sorkin 
conjectured that the expected value of its minimum /c-assignment is 

V 1 



' (m — i)(n — j) 

i,j>0, i+j<k v ' y J ' 

and they (with Aim) have proven this for k < 4 and in certain cases 
when k = 5 or k = 6. They were motivated by the special case of 
k = m = n, where the expected value was conjectured by Parisi 
to be Yli=i p- I n this paper we describe our efforts to prove the 
Coppersmith-Sorkin conjecture. We give evidence for the following 
stronger conjecture, which generalizes theirs. 

Conjecture Suppose that n, . . . ,r m and ci, . . . , c n are positive 
real numbers. Let X be a random m x n matrix in which entry Xij 
is chosen independently from the exponential distribution with mean 
— — . Then the expected value of the minimum ^-assignment of X is 

m + n - 1 - \I\ - \ J\\ 1 



n fc-i-|/|-|j| fm + n-l-\I\-\J\\ 
\ k-l-\I\-\J\ J 



Here the sum is over proper subsets I of {1, ... , m) and J of {1, . . . , n} 
whose cardinalities |/| and \J\ satisfy [/[ + \ J\ < k. 



* All three authors are affiliated with the Center for Communications Research, Prince- 
ton, NJ 08540. Email: buck@idaccr.org, clara@idaccr.org, robbins@idaccr.org 
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1 Problem Description and Background 



Suppose that k, m and n are positive integers with k < m < n. A minimum 
k-assignment of an m x n matrix X is a set of k entries of X, no two of 
which belong to the same row or column, whose sum is as small as possible. 
We denote the value of this minimum sum by mhifc(X). 

We say that a random real number x is exponentially distributed with 
rate a if it is chosen according to the density ae~ ax , x > 0. The mean value 
of a rate a exponentially distributed quantity is 1/a. 

Suppose that we generate a random m x n matrix X by choosing each 



entry independently from the exponential distribution with rate 1. In [ 05 
Coppersmith and Sorkin conjectured that the expected value of its minimum 
fc-assignment is 

Conjecture 1 

E(mm k (X))= Yl 1 r C 1 ) 

z — ' (m — i)[n — 7 

i,j>0, i+j<k v /v J ' 

In [|AS|| Aim and Sorkin show that this conjecture is correct when k < 4, 
when k = m = 5, and when k = m = n = 6. 

The conjecture of Coppersmith and Sorkin generalized a conjecture of 



Parisi [0] who considered the case k = m = n. In this CcLSCj clS shown in IICSI 
([I]) reduces to 

£(min fc (X)) = \r (2) 
i=i 

In this paper we describe our efforts to prove these conjectures. Our main 
result is Conjecture which generalizes the Coppersmith-Sorkin conjecture. 

We will say that a matrix X is random exponential with rate matrix A = 
{(lij) if each entry Xij is chosen independently according to the exponential 
distribution with rate a^ . The expected value of the minimum /c-assignment 
of such a matrix X is then a function of the rate matrix A. We denote this 
function by Ef.(A). 

We will show that Ek(A) is a rational function of the rates a^- and give an 
explicit method for computing it, at least in principle. Then we will specialize 
to the case when the rate matrix has rank 1, for which we have the following 
explicit formula. 
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Conjecture 2 Suppose that r 1; . . . , r m and Ci, . . . , c n are positive real num- 
bers and that = r^Cj . Let X be a random m xn matrix in which entry Xij 
is chosen independently from the exponential distribution with rate a^. Then 
the expected value of the minimum k-assignment of X is 

Vf-1 \k-i-\i\-\J\ (m + n-l-\I\-\J\\ 1 

Here the sum is over proper subsets I of{l, . . . , m} and J of {1, ... ,n} whose 
cardinalities \I\ and \J\ satisfy \I\ + | J\ < k. 

Example 1 The expected value of the minimum 1-assignment of a random 
exponential matrix with rate matrix 

l 



Example 2 The expected value of the minimum 2-assignment of a 3 x 3 
random exponential matrix with rate matrix 

1 1 1 
+ + 



^2 + ^3 n + r 3 ri+r 2 J c x + c 2 + c 3 

1 1 1 \ 1 

+ I : + : + 



c 2 + c 3 ci + c 3 ci + c 2 / n + r 2 + r 3 
5 



(ri+r 2 + r 3 )(ci + c 2 + c 3 ) 

We will provide evidence in support of Conjecture |2|. 

We also have a stronger conjecture for which we will provide evidence, 
although perhaps this evidence is not as strong as that for Conjecture [| A 
matrix can have several minimum ^-assignments for some value of k. How- 
ever, with probability 1, a random matrix has a single minimum ^-assignment 
for each k. Suppose that M is a (k — 1) x (k — 1) submatrix of X and that 
Xm{X) is the function with value 1 when M contains a minimum (k — 1)- 
assignment of X and otherwise. Then we define the expected contribution 
of M to the minimum /c-assignment of X as the expected value of the ran- 
dom variable XM{X)mmk(X). It is clear that E(mmk(X)) is the sum of the 
expected contributions of all the (k — 1) x (k — 1) submatrices. Our stronger 
conjecture gives a formula for the expected contribution of M when the rate 
matrix has rank 1. 



3 



Conjecture 3 Suppose that A = (rjCj) is a positive mxn matrix with rank 1 
and that X is a random exponential matrix with rate matrix A. Let I be a set 
ofk — 1 elements of {1, . . . , m}, let J be a set ofk — 1 elements of {1, . . . , n}, 
and let M be the (k — 1) x (A; — 1) submatrix of X with rows indexed by I 
and columns indexed by J. Then the expected value of xm(X) ■ min&pf) is 

Efio L ^Jtt ^t-t ) e 



where the outer sum is over permutations (ii, . . . , ik-i) and (ji, . . . ,jk-i) of 
I and J, respectively, and R and C denote the sums of all 's and all Cj 's, 
respectively. 

We shall see that Conjecture § implies Conjecture || and that Conjecture § 
in turn implies Conjecture |T]. 

Section discusses what we know for the expected minimum assignment 
when the rate matrix is arbitrary. 

In Section ^| we discuss the way we arrived at Conjecture || and give some 
equivalent formulations, one of which is directly implied by Conjecture 

We discuss the computational evidence for our conjectures in Section f|. 

Section || gives additional evidence for Conjecture §. 

We would like to thank Jim Propp for bringing this problem to our at- 
tention. 



2 Theory for a general rate matrix 

2.1 Expected value for a general rate matrix 

We begin by showing that the general formula for the expected value of the 
minimum assignment of a random exponential matrix is a rational function 
of the rates, with denominators factoring into linear terms of special form. 

Recall that k, m, n are positive integers with k < m < n and that A = 
i^ciij^j is a positive mxn matrix. We form a random matrix X by choosing 
independently from the exponential distribution with rate Oy. The expected 
value of the minimum fc-assignment of X is then a function of A, which we 
will denote by Ek(A). 
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By definition of expected value, Ek(A) is given by the integral expression 




Ek(A) = ( [ | ) / mm k (X)e- AX dX 
<x>o 

where the integral is taken over the space of all nonnegative matrices X. 
Here A ■ X denotes the dot product ^ . aijXij and dX denotes the product 

We denote by Si- the set of all m x n matrices a such that all the entries 
of a are O's except for k entries which are l's, no two in the same row or 
column. There are k\ (™) (?) such matrices in Sk and these we identify in the 
obvious way with the possible locations of the minimum k- assignment of X. 
In particular, 

mhifc(X) = min (a ■ X) . 

For each a we denote by P a the set of nonnegative matrices X for which 
minfc(X) = er ■ X\ that is, P a is the set of nonnegative matrices X for which 
the minimum /c-assignment is a. Thus, we have 



E k (A) 




[ (a ■ X)e~ A - x dX 



(3) 



Note that each of the sets P a is a polyhedral cone determined by a fi- 
nite set of homogeneous linear inequalities a ■ X < r • X for all r G Sk- 
As a consequence, each P a can be decomposed into a finite collection C a of 
simplicial cones. It seems difficult to give an explicit description of C a . Nev- 
ertheless, we can derive some useful properties of Ek(A) from the fact that 
this decomposition exists. First we rewrite ([3]) as 

Ek(A) = (n%) E E \f^-X)e- AX dX 



i,j / o-e<s fc ceC 



Each cone C is the set of nonnegative linear combinations of a set of mn 
linearly independent vectors Vi, i = 1, . . . , mn, where each V{ is a nonnegative 
m x n matrix. For the part of the integral over C, we make the substitution 
X = Ylii u iVii "where U = (ui, . . . , u mn ) ranges over all nonnegative mn- 
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tuples. We can then explicitly compute the integral over C as 



/ 1? ■ xy 

Jc 



~ Ax dX 



Ju>o i=1 

mn 

detV\J2 



mn 



U>0 



detV\ 



,i=i 



o-Vi 
A ■ Vi 



n 

vi=l 



where | det V| is the mn- volume of the parallelepiped determined by V±, . . . , V v 
Thus, we obtain the expression 

^-fe-Jss'^lS^)®^)- (4) 

Note that although the vectors Vi depend on C and a, they do not depend 
on A. Thus, Ek{A) is a rational function of the a^-'s, homogeneous of degree 
-1. 

We can obtain more information about the rational function Ek(A) by 
constructing, for each a G <Sfe, a finite set of generators of P CT in the sense that 
every element of P a is a nonnegative linear combination of the generators. 

For this purpose we define two classes of matrices. First, for any 1 < 
i < m and 1 < j < n we define to be the matrix that is all zero except 
for a single 1 at position Next, for any sets I C {l,...,m} and 

J C {1, . . . ,n}, we define Vj-j to be the matrix obtained from the all l's 
matrix by zeroing out all entries in the rows indexed by / and the columns 
indexed by J. It is easy to see that min^V/j) = max(0, k — \I\ — \ J\). Thus 
Vu is in P a if and only if a ■ Vu = max(0, k — \I\ — \J\). 

Theorem 1 Every element of P a is a nonnegative linear combination of 's 
with 6ij ■ o = and Vu 7 s in P a with \I\ + \J\ < k. 

We prove Theorem |I] using a reduction procedure on the matrices of P a . 
Let X be a matrix in P a and suppose that minfc(X) = s. We choose an arbi- 
trary linear order for the e^'s and denote this ordered set by ei, e 2 , • • • , e mn . 
Then we choose a sequence of nonnegative real numbers ax, ot 2) . . . , a mn as 
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follows. Once aci, . . . , are chosen we select aj as large as possible so that 
X — («iei + a 2 e 2 + • ■ • + ccjej) is nonnegative and has minimum fc-assignment 
with value s. 
Set 

Y = X - (aid + a 2 e 2 H h a mn e mri ). 

Note that if • ex 7^ 0, then we will have «j = 0, since otherwise a ■ 
(X — (ajei + a 2 e 2 + • — h a^e*)) < s. Thus X is K plus a nonnegative linear 
combination of the e^-'s given in Theorem |l|. 

We say that an entry of a matrix K participates in a minimum fc- 
assignment if there is a minimum /c-assignment using the entry y^. 

We say that a nonnegative matrix Y = (yij) is k-reduced if every nonzero 
entry of Y participates in a minimum ^-assignment. 

It is straightforward to see that the matrix Y resulting from our reduction 
process applied to X e P a is fc-reduced and that Y G P a . It remains to 
show that every &;-reduced matrix Y with minimum /c-assignment a is a 
nonnegative linear combination of the appropriate V}j. This will require a 
series of preliminary results. 

First we need a simple combinatorial lemma. 

Lemma 1 Suppose that T is a matrix all of whose entries are 0, 1, or 2 
and whose row and column sums are at most 2, and that the sum of all the 
entries in T is 2k. Then T = a + t for some k-assignments a and r. 

Proof: We may assume there are no 2's in T, since if there is a 2 we know that 
both a and r must have a 1 there, and only 0's everywhere else in its row and 
column. So we assume T is a 0-1 matrix whose row and column sums are at 
most 2, such that the sum of all entries is 2s for some s < k, and we want 
to find two s-assignments a and r such that a + r = T. Identify T with a 
graph with vertices at each 1 of T, and edges between any two l's belonging 
to the same row or column. Clearly every vertex of T has degree < 2, so 
every component of T is a chain or a cycle. The vertices in each component 
can be alternately assigned to o and r. If there is an odd component (which 
must be a chain) there must be another odd component to balance it out (so 
one can have an extra o vertex, and the other can have an extra r vertex), 
since the sum of all entries in T is even. □ 



7 



Lemma 2 Suppose that Y is a k -reduced matrix and 



S 



a b 
c d 



is a submatrix ofY. If a and d each participate in a minimum k-assignment 
of Y , then a + d < b + c. If also a + d = b + c, then b and c also each 
participate in a minimum k-assignment. These statements also hold for a 
and d switched with b and c. 

Proof: Suppose that a + d > b + c. Let (j\ and T\ be minimum assignments 
passing through a and d respectively. Form the matrix T\ = o~\ + T\. Then 
form T from T\ by subtracting 1 at the positions of a and d and adding 1 
at the positions of b and c. The hypotheses of Lemma [I] apply to T so that 
T = a + t for some fc-assignments a and r. But since a + d>b + cwe must 
have 



contradicting the minimality of the assignments o\ and T\. Thus a + d < b+c. 
Now suppose that b + c = a + d. Then the same construction yields 



so that both a and r are minimum /c-assignments, and at least one includes 
b and at least one includes c. 

This proof obviously also holds with a and d switched with b and c. □ 

Proposition 1 Suppose that Y = (yy) is a k-reduced m x n matrix. Then 
there exist Ax, ... , A m and fii, . . . , \x n such that 



and such that yij participates in a minimum k-assignment precisely when 
Vij l>, ■ Aj. 

Proof: Let d = yt u denote the largest entry in Y . Take Aj to be the i th entry 
in the column of d, so that Aj = y iu . Let fij to be the j th entry in the row of 
d, decreased by d, so \ij = y t j — d. 
First we prove @. 



a 1 -Y + r 1 -Y = T 1 -Y>T-Y = o--Y + T-Y, 



a 1 -Y + r 1 -Y = T 1 -Y = T-Y = a-Y + T-Y, 



y i:j = max(0, Aj + fij) 



(5) 
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When yij is in the row or column of d, then we have y^ = Aj + fij, so (^) 
is immediate. 

Suppose that a is in neither the row of d nor the column of d. Let 

S — a ^ 
c d 

be the submatrix of Y containing the rows and columns of a and d (where 
the order of the rows or columns in S may be opposite to the order they 
occur in Y). Then we just need to show that a = max(0, b + c — d). 

First suppose that b + c — d > 0. Then b, c, d must all be positive, because 
d is maximum for the whole matrix. Since b, c are both positive, they must 
both participate in a minimum k- assignment, so b + c < a + d, by Lemma |2|. 
But then, a > b + c — d, so a is positive. Then a and d both participate in a 
minimum k- assignment, which implies a + d < b + c, so a = b + c — d. 

Next, suppose that b + c — d < 0. Ifa>0, then a and d both participate 
in a minimum /c-assignment, so d<a + d<b + c, a contradiction. Thus 
a = 0. 

Now we show that = Xi+fij exactly when y^ participates in a minimum 
k- assignment. 

Recall that for any y^ in either the row or column of d, we have y^ = 
Aj + fij. So we need to show that all such entries participate in a minimum 
k- assignment. Let b be any entry in the column of d. We already know that 
positive entries must participate in a minimum /c-assignment, so we assume 
that 6 = 0. If d = 0, then our whole matrix is zero and our result is trivial, 
so we may assume that d > and therefore participates in a minimum k- 
assignment. If the minimum assignment using d does not use the row of b, 
we can replace d by b and obtain a smaller assignment, a contradiction. So 
we can conclude that the minimum assignment using d also uses an element 
a from the row of b. Form the 2x2 submatrix S containing a and d as above. 
Since d is the largest entry, c < d, and we also have b = < a. Thus we can 
exchange a and d for b and c, to obtain a minimum fc-assignment in which 
b participates. In the same way, we see that any entry in the row of d must 
participate in a minimum /c-assignment. 

Finally, consider an element a that is neither in the row nor the column 
of d and form the 2x2 submatrix S containing a and d as above. We must 
show that a participates in a minimum assignment exactly when a = b+c—d. 
This is certainly true if a > 0. So, let us assume that a = 0. Also, since 
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both b and c are in a row or a column of d, both participate in a minimum 
assignment, so that b + c < a + d. 

Now suppose that a participates in a minimum k- assignment. Then a + 
d < b + c, so a = b + c — d, as required. 

Conversely, suppose that a = b + c — d. Then, since b and c participate 
in minimum k- assignments, Lemma |2| shows that a also participates. □ 

Proof of Theorem [J: Now we are ready to prove Theorem [l] by showing 
that any fc-reduced matrix Y in P a is a nonnegative linear combination of 
a suitable collection of matrices Vu from P a . Without loss of generality we 
can assume that the A's and /i's are weakly increasing. In this case the rows 
and columns of Y are also weakly increasing. 

If the matrix Y is zero, there is nothing to prove. Otherwise, since every 
nonzero entry in Y participates in a minimum fc-assignment , we know that 
the minimum fc-assignment is nonzero. Hence there is at least one nonzero 
entry among yi t k, y2,k-i, • • • , Vk,i- 111 particular, there is a pair such that 
Uij > and i + j < k + 1. Now select such a pair to be minimal in 
the sense that if % ^ 1 then Vi-ij = and if j ^ 1 then y%^-\ = 0. Let 
I = {1, . . . , i — 1} and let J = {1, . . . , j — 1}. We will show that Vu is in P a 
and that Y — yijVu is again in P a and still /c-reduced. 

Suppose that 1 < %' < i and 1 < j' < j. Then, since 

max(0, Aj/ + jUj) = yi'j = < y^ = max(0, Aj + yUj), 

we know that Aj/ < Aj. Similarly \iy < \iy Since y^j = max(0, A^ + fij) = 0, 
we have Aj/ + fij < 0. But then 

Aj/ + /ij/ < Aj/ + yUj < 0. 

If follows from Proposition |l] that none of the matrix entries y^y with i' < i 
and j' < j can participate in a minimum fc- assignment. 

To see that Vjj G P CT , first note that any minimum fc-assignment of Y 
must use all of the first i — 1 rows. If not, since % < k, there is some %\ > % 
such that row i\ participates. Then, since not all of the first i — 1 rows are 
used, we can replace the entry of the assignment in row %\ with the entry 
in the same column of row z , for some z < i, to get an assignment with 
a value no larger. The entry being replaced could not come from a column 
preceding j, by the discussion above, so it must be positive. But also by 
the discussion above, Aj < A^. Then the new assignment would be strictly 
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smaller, a contradiction. Thus, any minimum assignment uses all of the first 
i — 1 rows and all of the first j — 1 columns, and does not use any entry 
which is in both the first i — 1 rows and the first j — 1 columns. So, if 
r is any matrix representing a minimum fc-assignment of Y, we must have 
Vu ■ t = k — \I\ — \ J\. In particular, Vu ■ a = k — \I\ — \J\, so Vjj is in P a 
as claimed. 

The preceding argument shows that if we replace Y by Y' = Y — tVu, 
for any t satisfying < t < yij, the effect on any minimum fc-assignment is 
to subtract (k — \I\ — \J\)t from its value. Thus, all assignments r that are 
minimum for Y will agree on Y' . We now show that each of these assignments 
r is minimum for Y' as well. Assume not. Then we could find t\ and t 2 such 
that < ti < t 2 < Vij and A;- assignments r and such that r is minimum 
for (Y — £Vjj) when t < t\ but not when t\ < t < t%, and <fi is minimum for 

Y — tVu when t\ < t < t^. Thus, we would have both <fi and r minimum 
assignments for Y — t\Vu. But then our preceding argument applied to 

Y — tiVjj tells us that <fi and r must agree on Y — tVu for £i < t < t 2 , a 
contradiction. 

Now we let Y' = Y — yijVjj and observe that Y' is fc-reduced. Indeed, all 
the assignments r that were minimum for Y are also minimum for Y' . Thus, 
any element that participated in a minimum /^-assignment for Y will also 
participate in a minimum /c-assignment for Y' . Also, the replacement of Y 
by Y' creates no new nonzero elements, so the matrix Y' will be A;-reduced. 

Since a in particular is a minimum assignment for Y, a will also be 
minimum for Y', so that Y' is also in P a . 

Note that Y' has at least one more zero entry than Y, namely the entry 

at 

We can continue removing multiples of submatrices Vu, each time pro- 
ducing a matrix with at least one more zero entry. Thus, we eventually reach 
a matrix that is all zero. In effect, we have expressed Y as a nonnegative 
linear combination of the generators as required. □ 

We have shown that every element of P a is a nonnegative linear com- 
bination of certain e^'s and Vjj's in P a . We remark that the e^'s generate 
extreme rays of P a but the Vi/s in general do not. It is not hard to show that 
the Vijs which do generate extreme rays of -Pu are those with |J| + |«7| = k— 1. 

We now observe that Theorem [I] allows us to make some conclusions 
about the rational function Ek(A). A simplicial cone in a decomposition of 
P a has some generators of the form and some of the form Vjj. For the 
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generators of the form e^, we know that a ■ ey = 0. Also the dot products 
eij ■ V in the denominator cancel with factors in the initial product of o^'s in 
d^). Thus, the denominator of the integral over a simplicial cone is a product 
of terms of the form A ■ Vij. Finally, we can conclude: 

Theorem 2 The expected minimum k-assignment of a random exponential 
matrix with rate matrix A is a rational function of the entries of A. The 
denominator of the rational function is a product of sums, each being the 
sum of all entries in a submatrix of A omitting i rows and j columns, where 
i + j < k. 

Example 3 The minimum 2-assignment of a random 2x2 exponential ma- 
trix with rate matrix A = (a^-) is 

1 0-11^21 

a 2 i + a 22 (an + ai 2 )(a 2 i + a 22 )(ai 2 + a 22 ) 
Gl 2 Q 22 

(an + a 12 )(a 21 + a 22 )(an + a 2 i) 

The above formula is easily computed by the method we will sketch in Sec- 
tion f§ 

In the case of a rank 1 rate matrix A = (rjCj), Theorem |2| says that 
Ek(A) is a rational function of the rVs and c/s and its denominator is a 
product of sums of subsets of the ?Ys omitting fewer than k of the r,'s and 
sums of subsets of the c/s omitting fewer than k of the Cj's. Since this is a 
consequence of Conjecture 0, it lends some support to the conjecture. 

2.2 The nesting lemma 

A real nonnegative m x n matrix X has minimum fc-assignments for each 
k < m. Generically there is only one of each but in some cases there are 
many minimum assignments of various sizes. It helps to know how these are 
related. 

The following lemma is fundamental. Other proofs probably exist but we 
include ours here for completeness. 

Lemma 3 Let k\ and k 2 be two integers, with ki < k 2 < m. 



1 



All + Oj\ 2 



+ 



+ 
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Suppose that M 1 is a k\ x k\ submatrix of X that contains a minimum 
ki-assignment of X. Then there exists a k 2 x k 2 submatrix M 2 containing 
Mi such that M 2 contains a minimum k 2 -assignment of X . 

Suppose that M 2 is a k 2 x k 2 submatrix of X that contains a minimum 
k 2 -assignment of X. Then there exists a k± x k\ submatrix Mi contained in 
M 2 such that M x contains a minimum ki-assignment of X . 

Proof: If ki = k 2 , there is nothing to prove. So assume ki < k 2 and fix a 
minimum /ci-assignment and a minimum /^-assignment. 

Let G be the graph on k\ red vertices (representing the entries of the ki- 
assignment) and k 2 blue vertices (representing the entries of the /^-assignment), 
with edges between two vertices if the corresponding entries belong to the 
same row or column. (If the assignments share an entry then we have a 
red vertex and a blue vertex with two edges between them comprising a 
component which is a cycle of length 2.) 

Then G is bipartite and no vertex of G has degree more than 2. Thus, 
every component of G is a cycle or a chain in which the red and blue vertices 
alternate. 

Suppose some component of G has mi red vertices and m 2 blue vertices, 
and all of its red vertices have degree 2. Such a component is either a cycle, so 
that mi = m 2 , or a chain with blue vertices at each end, so that mi + 1 = m 2 . 
In either case we have m\ + 1 > m 2 . 

Consider the submatrix Mi spanned by the associated mi entries of the 
/ci-assignment and the submatrix M 2 spanned by the associated m 2 entries 
of the ^-assignment. The component condition translates into the condi- 
tion that the Mi is contained in M 2 . Also the remaining entries of the 
two assignments comprise a minimum (ki — mi)-assignment and a minimum 
(k 2 — m 2 )-assignment of the submatrix of X complementary to M 2 . Since 
ki < k 2 and m x + 1 > m 2 , we have k x — m x <k 2 — m 2 , so the lemma follows 
by induction applied to the complementary submatrix. 

Now assume that every component of G has a red vertex of degree 1 or 
less. Such components are chains with one endpoint red. When the other 
endpoint is red, there are more red than blue vertices in the component. 
When the other endpoint is blue, the number of vertices of both colors is 
equal. In particular the number of red vertices is always at least as great 
as the number of blue vertices. It follows that, all together, there are at 
most ki blue vertices that are connected to some red vertex. Thus we can 
select k 2 — ki entries of the /c 2 -assignment that do not share any row or 
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column with the /ci-assignment . We can now consider three sets of matrix 
entries, the set Si of ki entries of the /^-assignment, the set S 2 , just selected, 
of k 2 — k\ of entries of the ^-assignment and the set S3 of the remaining 
k\ entries of the k 2 - assignment. Then the sets S3 and Si are both in the 
submatrix complementary to that determined by S 2 . Moreover the sums of 
the entries in Si and S3 must be equal. For, if the sum of the Si entries were 
greater than that of S 3 , then S 3 would be a smaller /^-assignment than Si, a 
contradiction. But if the sum of the Si entries were smaller than the sum of 
S3 entries, then Si U S 2 would be give a smaller /c 2 -assignment than S 2 U S3, 
which was a minimum /^-assignment. 

Thus S3 is a minimum /^-assignment contained in our minimum k 2 - 
assignment and Si U S 2 is a minimum /^-assignment containing our minimum 
fci-assignment, which proves our lemma. □ 

2.3 A computational consequence of the nesting lemma 

Lemma |3| suggests a way to compute the minimum /c-assignment of a matrix 
X when k < n. We proceed by finding the minimum /c-assignment s for X 
for k — 1, 2, ... , one at a time. Suppose we have found a minimum (k — 1)- 
assignment whose rows and columns determine a (/c— 1) x (A;— 1) submatrix M 
of X. Then, when we search for a minimum k- assignment, we know that we 
can restrict our search to minimum /^-assignments in one of the submatrices 
of X obtained by appending a single new row and new column to M. 

There are some simple properties that such an extension must have. Sup- 
pose the submatrix of a minimum /c-assignment uses M together with a new 
row i and a new column j. Also, suppose that the minimum /c-assignment 
uses an entry from column j that is in row %' of M. Then this entry in column 
j must be the smallest entry in the part of row i' outside of M. Similarly, if 
the minimum /c-assignment uses an entry from row % that is in column j' of 
M, then this entry in row % must be the smallest entry in the part of column 
j' outside M. Finally, if the minimum /c-assignment uses the entry Xy, then 
this entry must be minimum in the submatrix of X complementary to M. 

The preceding discussion shows that the following strategy will work to 
construct the minimum /c-assignment once we have found the minimum (k — 
l)-assignment. We define the kx k auxiliary matrix Auxm{X) by appending 
to M a new row and column as follows. To each row i of M we append a 
new entry which is the minimum of the entries in row i of X that are outside 
of M. To each column j of M we append a new entry which is the minimum 
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of the entries in column j of X that are outside of M. At the intersection of 
the new row and new column we place the minimum entry of the submatrix 
of X complementary to M. 

We now find the minimum /c-assignment of A\xxm{X). This assignment 
will use an entry in the last row and an entry in the last column of Auxm(A), 
which can be the same entry if the assignment uses the entry in the last row 
and column. Each of these entries is a copy of some entry of X, which then 
tells us which row and column of X need to be appended to M to obtain the 
submatrix of the minimum /c-assignment of X. 

2.4 A consequence for the expected contribution 

The discussion of the auxiliary matrix in the preceding section can be for- 
malized to prove an interesting property of the expected contribution. 

We will use the following simple facts about independent exponential 
random variables. 

Proposition 2 Suppose that a±, . . . , a m are positive real numbers, and that 
Xi, . . . ,x m are independent random variables with Xi chosen from the expo- 
nential distribution with rate Oj. 

Let x denote the random variable minjXj. Then x is distributed as an 
exponential random variable of rate a\ + ■ ■ ■ + a m . 

Let W be the discrete random variable whose value is the least i for which 
x = Xi. Then W and x are independent random variables, and the probability 
that W = i is a^j [a\ + • • • + a m ) . 

Proof: Let c be a positive real number. Then the probability that x > c is 



where we have made the substitution Xi = u, L + c on the second line. Taking 
the derivative with respect to c we see that x is distributed as an exponential 
random variable of rate a\ + • • • + a m . 
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A similar substitution yields the part of the integral corresponding the 
event that X — X ^ - 3jS follows: 



= a\ ■ • ■ a, 



a\ ■ ■ - a 



III 



111 




= 



e 




(ai (c+ui )+a 2 (c+«i +u 2 )+- ■ -a m (c+«i +u. 



e 



1 " ' ■ c/./'| • • • dx, 



>^dui ■ ■ ■ du. 



■in 



Em 



e 



Thus, the conditional probability that x = Xi, given x > c, is ai/{a\ + • • • + 
a m ). Since this is true for all c, the event X — x^ IS independent of the random 
variable x. □ 

Now let X be any nonnegative mxn matrix and associate to X the k x k 
matrix Y = Auxm(AT) where M is the upper left (k — 1) x {k — 1) submatrix 
of X. For any t < k, let M t denote the upper left t x t submatrix of X. By 
abuse of notation we will also let M t denote the upper left t x t submatrix 
of Y, since they are identical. 

Lemma 4 For any t < k, let a be a t- assignment of M t . Then a is a 
minimum t-assignment for X if and only if it is a minimum t-assignment 
for Y . 

Proof: Let a be a t-assignment of M t . 

If t — 1 it is clear that the lemma holds. We proceed by induction on t. 

Suppose that a is a minimum t-assignment of X. By induction on t, there 
is a minimum (t — l)-assignment of Y which lies in M t , so by Lemma there 
is some minimum t-assignment r of Y which uses at most one column and 
one row outside of M t . By definition, each entry of Y outside of M t equals 
some entry of X outside of M t . If we replace each of the entries of r outside 
of M t by their equivalent entries in X, then we get a t-assignment r' of X. 
Furthermore, r' • X > a ■ X, since a is minimum. But r' ■ X = r • Y and 
a ■ X = a ■ Y , so r ■ Y > a ■ Y . Therefore, a is a minimum t-assignment of Y. 

Now suppose that a is a minimum t-assignment of Y. By induction on t, 
there is a minimum (t — l)-assignment of X which lies in M t , so by Lemma |3|, 
there is some minimum t-assignment r of X which uses at most one column 
and one row outside of M t . We can replace each entry Xij of r by the entry 
yioia(i,k),mm(j,k) to get a t-assignment r' of Y. Furthermore, r' • Y > a ■ Y, 
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since a is minimum. But r • X > r' ■ Y and a ■ X = a ■ Y, so r • X > a ■ X. 
Therefore, a is a minimum t-assignment of X. □ 



Lemma 5 The minimum (k — l)-assignment of X uses its first k — 1 rows 
and columns exactly when the minimum (k — 1) -assignment of'Y uses its first 
k — 1 rows and columns. In this case, the value of the minimum k-assignment 
of X is the same as the value of the minimum k-assignment ofY. 

Moreover in this case there is a minimum k-assignment t of X and a 
minimum k-assignment r' of Y which correspond entry by entry; i.e., the 
entries of t in M are located in the same positions as the entries of r' in 
M , and for any entry Xij of r outside of M there is a corresponding entry 

Vinm(i,k),mm(j,k) of V Outside of M . 

Proof: The first paragraph is the case t = k — 1 in Lemma ^. The second 
paragraph follows from our discussion of auxiliary matrix in the preceding 
section. □ 

We will use the more detailed statement in the second paragraph to prove 
Lemma ^| in Section |372. 



Now let X be a random exponential m x n matrix with rate matrix A, 
and associate to A the k x k matrix B = (bij) with 

bij = dij when 1 < i, j < k — 1 

and 

kk = ^2 ai i' ' i = 1, ■ ■ ■ ,k - 1 

j>>k 

and 



i'>k 



and 



%k — a i'j'- 
i',j'>k 



Note that Y = Auxm(X) is a random exponential matrix with rate matrix 
B, by Proposition |[ 

The statement in the first paragraph of the preceding Lemma shows that 
the expected contribution of the submatrix of X consisting of its first k — 1 
rows and first k — 1 columns to the minimum /c-assignment of X is the same 
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as the expected contribution of the submatrix of Y consisting of its first k — 1 
rows and first k — 1 columns to the minimum /^-assignment of Y. 

In the rank 1 rate matrix case, if C(k, 7*1, ... , r m , c\, . . . , c n ) denotes the 
expected contribution of the submatrix of the first k — 1 rows and columns 
when the rate matrix is A = (riCj), then 

C(k, ri, . . . , rfc_i, r k H h r m , Ci, . . . , c k -i, c k H h c n ) 

is the expected contribution when the rate matrix is B. These two functions 
must agree, so the entries of A outside of the first k — 1 rows and columns 
enter into the contribution function only via the sums defining the entries in 
the k th row and column of B. 

Now note that Conjecture has the feature that it predicts the same 
contribution in the two cases above. So this is a bit of evidence in favor of 
Conjecture 0. 

Thus, we can summarize our discussion in the rank 1 rate matrix case as 
follows: 

Theorem 3 Let k, m and n be integers k < m < n. Then Conjecture [| 
holds for k- assignments in an m x n matrix if and only if it holds for k- 
assignments in a k x k matrix. 

Thus, if one could prove only the cases k = m = n of Conjecture [| that 
would prove the general case of that conjecture as well as Conjecture ^] and 
Conjecture [I]. 

3 Rank 1 rate matrices 

In this section we restrict our attention to random exponential matrices for 
which the rate matrix has rank 1. We describe how we arrived at Conjec- 
ture and then give equivalent formulations, which provide different kinds 
of confirmation. Finally, we prove a probability result about the locations of 
the minimum ^-assignments for 1 < £ < k. 

We discovered Conjecture |2| while experimenting with the computations 
described in Section f|. Finding that Mathematica had trouble carrying out 
the computation when the rate matrix consisted of mn indeterminates, we 
decided to try a simpler case, with rate matrix of the form = aj for all i and 
j. We noticed in this case that the answer has a surprisingly simple form — in 
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particular, it can be written as a linear combination of the reciprocals of 
sums of the a,'s. Next we found that, when the rate matrix has rank 1, so 
that ciij = TiCj, the expected value seems to be a linear combination of terms 
of the form 

(6) 



with / a proper subset of {1, . . . , m}, and J a proper subset of {1, . . . , n}. 
It is easily shown that the rational functions 7^ ^ r are linearly 

independent over the real numbers. Thus the coefficients in such a linear 
combination are uniquely determined. 

Making the assumption that the expected value is indeed a linear com- 
bination of terms of the form (|6]), we arrived at Conjecture |2| by considering 
certain limiting conditions on the expected value. We will describe these 
limiting conditions in Section |5|. 

3.1 Equivalent formulations of Conjecture [2] 

From now on we use the shorthand notation [m] for the set {1, . . . , m}. 
Let us introduce the notation 

m + n — 1 — III — I J\\ 1 



^,^B- 1 )™-( m r\\ i r/ i r J r i ) 

i, j 



I\-\J\ >/(£*jr,) 

(7) 

for the formula given in Conjecture |^. Here recall that r = (n, . . . ,r m ) is 
an m-tuple of positive real numbers, c = (ci, . . . , c n ) is an ra-tuple of positive 
real numbers, and the sum is over proper subsets I C [m] and J C [n]. The 
binomial coefficient enforces the condition |/| + \J\ < k. In what follows we 
will often not mention such constraints explicitly. 

In this section we derive alternative ways to write (I7|) and conclude that 
F(k,r,c) is positive, Conjecture || implies Conjecture [l], and Conjecture § 
implies Conjecture 0. 

Note that F(k,r,c) can be written more succinctly as 



W>,r,c) = E( fc _!_ I/, _|j|J 



1 



A) (E^/ r i)(E^j^) : 



(8) 



using binomial coefficients with negative numerator. 
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Proposition 3 

F(t,r lC )= Y. (-l)^Mi.HH^H^ ___i_ _ . (9) 

|/'| + |J'|<fc,/C/',JCJ' 



i/ere J' and / are proper subsets of \m\ and J' and J are proper subsets of 
[n\. 

Proof: Comparison with (||) shows that, for a fixed / and J with |/| + |J| < k, 
we need to evaluate the sum 

(_i)l*'H-WI-kl. 

ICI',JCJ',\I'\+\J'\<k 

If we denote by i and j the cardinalities of I and J and by t and u the 
cardinalities of I' and J', we can rewrite this sum as 



E (-D 

<>j,ti>jr,t+U<fc 



t-i )\u- j 



t>0,n>0,i+u<fc-i-j 
fe— 1— i— j 

E (-1)' E 

i=0 t>0,u>0,t+u=l 
k — l—i—j 



m — i\ I n — j 
t ){ u 

m — i\ ( n — j 
t )[ u 



E~(-i) , ( w "~ i ~ J ) 



1=0 

k—l—i—j 

E 

1=0 

k—l—i—j 

E 



1=0 

k — 1 — m — n 
k—l—i—j 



which agrees with (H). □ 



i + j — m — n + I — 1 
I 



i + j — in — n + l\ fi + j— m — n + l — 1 
I ~\ l-l 
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We can rewrite as a double sum with the inner sum over / and J and 
the outer sum over I' and J'. Then, for fixed I' and J', the inner sum factors 

as 




Now we show that each factor has an interesting probabilistic interpretation. 

Suppose that an urn contains m balls labeled 1,2, ... ,m and for each 
i, ball i has weight r;. We select balls one at a time without replacement, 
at each time selecting a ball with probability proportional to the weights of 
those balls still in the urn. Let Pr(r, I') denote the probability that the set 
of balls in /' are the first t balls to be chosen, where t is the cardinality of I'. 
Then 

7T i = l ^ 2-Uj=l 

where R = YlT=i r « anc ^ the outer sum is over ai l ^ orderings (tti, . . . , ir t ) of 
/'. 

We can calculate Pr(r, I') in a different way as follows. Suppose we draw 
all m balls from the urn. If we fix any subset U of balls, then the probability 
that a particular ball u from U is chosen before any other ball from U is the 
weight of u divided by the sum of the weights of the balls in U. 

Now, for i G let Ei denote the event that the first time a ball is drawn 
from the set consisting of i together with the complement of I', the ball 
chosen is from the complement of Then Ei has probability 

In order for our set I' to be the set of the first t balls chosen, it is necessary 
and sufficient that none of the events E^ occur. For any subset I of I' the 
probability that all of the events E iy i G /, occur is 
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So, by the Inclusion-Exclusion Principle, 



Pr(r, 



I^ = Yt-l)\r\-\i\Z&rIi (12) 



icr 



which is (Eigr r «) times the first factor in (|H]). The analogous result holds 
for the second factor in fllPp. We conclude that 

Proposition 4 

/,J,|/| + |J|<fc ^i^ r 'A^J C 3J 

where the sum is over proper subsets I of [m] and J of [n]. Hence F(k, r, c) 
is always positive. 

□ 

Now we rewrite F(k, r, c) to show that Conjecture ^| implies Conjecture p] 
and Conjecture |] implies Conjecture ||. 

Let Pr(r, (i 1; . . . , i t )) denote the probability that the first t selections from 
our urn are ii, . . . , i t in that order. Then, from (|l"3"D and (0), we get 

where the sum is over all sequences i = (ii, . . . , i t ) of distinct integers in [in] 
and j = (ji, . . . ,j u ) of distinct integers in [n]. But, certainly 

Pr(r, (ii, . . . , it)) = Pr(r, (i 1} . . . , i t , . . . , i fc _i)) 

where the sum is over all extensions of (ii, . . . , i t ) to a (k — l)-long sequence 
i of distinct integers in [m]. Thus, we can rewrite (|14]) and obtain: 

Proposition 5 

F(k,r, C ) = Y, E Pr(r,(2l ;"-'; fc - l))PrMjl :--- ,Jfc - l)) , (15) 



where the outer sum in ( fiq j zs over pairs of ordered sequences ofk — 1 distinct 
integers from [m] and [n] . 
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Note that each term in the above sum corresponds to a flag of submatrices 
of sizes 1 xl, . . . , (k — 1) x (k — 1). In this form, specializing to the case that all 
the r's and c's are 1, it is easy to see that Conjecture implies Conjecture [I]. 

We now group the terms in the outer sum according to the (unordered) 
sets / = . . . , ik-i} and J — . . . ,jk~i}- It then becomes 



F(k,r,c) = J2i: E ^^•••; tfc - lJJ ^u...,3k-i)) (16) 

I.J U t.u>0.t+u<k ( R Y,s=i r is)(C J2s=l C js) 



where the outer sum is over sets / and J of size k — 1 and the inner sum is 
over permutations . . . , ik-i) of I and permutations (ji, . . . ,jk-i) of J. In 
this form we can see that the term of the outer sum corresponding to the sets 
/ and J is the expected value of the contribution of the submatrix with row 
indices / and column indices J predicted by Conjecture |[ Since the sum 
of the expected contributions of all submatrices is the expected minimum 
/c-assignment, we now see that Conjecture |3| implies Conjecture [2|. 

Finally, for any T C I C [m], let Pr(r, T, I) denote the probability that 
the first |T| balls drawn from the urn comprise the set T and that the first 
|/| balls drawn comprise the set /. Then we can rewrite our formula for the 
expected contribution of a submatrix with rows / and columns J as 



E 



Pr(r,T,/)Pr(c,[/,J) 



TCI,UCJ,\T\ + \U\<k Etj£T r *)(E^t/ °u 



3.2 Flag probabilities 

In this section, we prove ap probability result in the special case that the 
rate matrix has rank 1. This result may be the reason that simple formulas 
exist for Ej.(A) when the matrix A has rank 1. 

It is possible for a matrix to have many minimum ^-assignments for some 
k. However, with probability 1, a random matrix X has a unique minimum 
k- assignment for each k. So, if we let denote a k x k submatrix of X 
containing a minimum fc-assignment of X, then, with probability 1, is 
unique. By Lemma ^, the submatrices are nested: Mi C Mi C • ■ ■ C 
Mfc. We will call this the flag of submatrices of X. This flag can also be 
described by the list ii, i 2 , ■ ■ ■ ,ik of appended rows and the list ji, j 2 , ■ ■ ■ ,jk 
of appended columns. Thus, Mi is the submatrix with rows and 
columns 
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It is natural to ask for the probability that a random matrix has a given 
flag of submatrices. We know of no formula for this probability for general 
rate matrices. However, for rate matrices of rank 1, we can prove a simple 
formula for the probability of each flag. Moreover, this formula will involve 
the probabilities Pr(r, (ii, . . . ,ik)) calculated in the previous section. 

We first need the following 

Lemma 6 Let X be an exponential random matrix with rank 1 rate matrix 
A = {riCj). If the minimum (k — 1)- assignment of X uses the first k — 1 rows 
and first k — 1 columns, then the minimum k-assignment uses an additional 
row and column. The probability of row i' being the additional row is — , 

the probability of column j' being the additional column is ^f 2 ^ — , and these 
events are independent. 

Proof: Let M denote the upper left (h — l)x(k — 1) submatrix of X, and 
Y the k x k matrix Auxm{X) defined in Section [2.3| . Then the upper left 
submatrix of Y is identical to M, and, by abuse of notation, we will let M 
denote that submatrix of Y also. If the minimum (k — l)-assignment of X 
lies in M, then by Lemma |], the minimum (k — l)-assignment of Y lies in 
M, and by Lemma [5], the minimum k- assignments of X and Y correspond 
entry by entry. 

There are two cases to consider. 

In the first case, for some s < k — 1 and t < k — 1, the minimum k- 
assignment of Y uses the entries y s k and yu- These entries correspond to 
the minimum entry in row s of X outside of M and the minimum entry in 
column t of X outside of M. In the second case, the minimum ^-assignment 
of Y uses entry y^. This entry corresponds to the minimum entry in the 
submatrix of X complementary to M. In both cases, by Proposition ^ the 
locations of these minima in X are independent of the random variables 
making up the entries of Y, and thus independent of the events that the 
minimum (k — l)-assignment of Y lies in M and the minimum /c-assignment 
of Y uses particular entries outside of M. 

Thus in the first case, the probability that the minimum entry in the part 
of row s outside of M comes from column j' is ^n sCj/ — = — r 3 -, — i and the 

l^j = k r = C 3 C fe + ... + C n ' 

probability that the minimum entry in the part of column t outside of M 
comes from row i' is ^li' Ct — = — ^ — . Moreover, the locations of these 

2^ i=k riC t r k + ...r m ' 

minima within row s and column t are independent events since the parts of 
row s and column t outside of M are disjoint. 
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In the second case, the probability of the minimum entry in the submatrix 
of X complementary to M coming from row %' and column j' is where 
the sum in the denominator is over all locations in the submatrix of X 
complementary to M. Thus, the probability that the minimum entry comes 
from row i' is — and the probability that the minimum entry comes from 

column f is — . □ 

From this lemma we can immediately conclude the following theorem, 
which imparts further meaning to the formal Pr(r, J) and Pr(r, (zi, . . . , ij.)) 
functions used in Section |3TT| : 

Theorem 4 Suppose that A = (rjCj) is a rank 1 rate matrix and X an 
exponential random matrix with rate matrix A. Let (ii, . . . , i^) be a sequence 
of distinct elements of [m] and (ji, ■ ■ ■ ,jk) a sequence of distinct elements 
from [n]. Then the probability that X has the associated flag of submatrices 
is Pr(r, (ii, . . . , i^)) Pr(c, . . . ,jk))- Furthermore, if I C [m] and J C [n] 
are sets of size k, then the probability that the minimum k-assignment of 
X uses the rows indexed by I is Pr(r, J) and the probability that it uses the 
columns indexed by J is Pr(c, J), and these events are independent. 

□ 

We will see that this formula enters in an essential way into the proof of 
Theorem 



4 Computational evidence for our conjectures 

By Proposition |2], it is an easy matter to compute the expected value of the 
minimum 1-assignment for an arbitrary rate matrix. 

Example 4 The expected value of the minimum 1-assignment when the rate 
matrix is A = (a^) is 

1 



For k > 2 the computation is more complicated. In [|AS|| and ||CS|| the au- 
thors calculate the expected value of the minimum assignment for a random 
exponential matrix when the rates are all 1 and k is small. 
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The method in | AS|| applies just as well to the case of arbitrary rate 



matrices. The essence of their idea is to introduce a slightly more general 
expectation problem in which they choose all the entries of the random matrix 
X as before, except that there is a set Z of fixed zeroes in X. Let us denote 
the expected value of the minimum assignment in this case by E(A, Z). 

It is then sometimes possible to establish a recursive calculation of E(A, Z). 
The base of the recursion occurs when there exist k zeroes in Z, no two in 
the same row or column. In this case we know that the expected value of the 
minimum assignment is zero. For the inductive part of the calculation we 
can sometimes express an expected value E(A, Z) as a constant plus a linear 
combination of E(A, Z') where Z' is obtained from Z by adjoining one more 
position to Z . 

This arises as follows. Suppose that X is a random exponential matrix 
except for a set Z of positions in X where the entries are fixed zeroes. Suppose 
further that we have a set S of positions in X, disjoint from Z, such that any 
minimum /c-assignment of X meets the set S in exactly r positions. (In other 
words, every nonnegative matrix with zero set Z has the property that its 
minimum fc-assignments all meet S in exactly r positions.) Abusing notation, 
we also let 5* denote the matrix which is 1 at the positions in the set S and 
zero otherwise. 

We will derive the following formula: 

E(A,Z) = -^+ ^E(AZU{(i,j)}). (18) 

Indeed, the integral for E(A,Z), which involves only the variables Xij for 
i Z, is given by 




E(A,Z)= | |7 dij ] / mm k (X)e~ A - x dX. 



We can derive ( 18 ) by breaking up this integral into \S\ parts, each corre- 
sponding to a position in S containing the minimum entry among all positions 
in S. For the part of the integral where Xi j is the minimum entry in S, we 
make a change of variables with Jacobian 1, as follows. We express the Xjj in 
terms of new variables yij by setting x^ = yi j +yij when G S — {(iq, jo)} 
and = otherwise. X can then be written as Y + yi j S where Y is a 



26 



nonnegative matrix with fixed zeroes at Z U {(io,jo)}- From our hypothesis 
about S, we have 

mm k (X) = min fc (F) + ry iojo . 

(Otherwise, there would be a non-minimum fc-assignment of X, meeting S 
in fewer than r positions, that becomes a minimum /c-assignment of a matrix 
X — tS for some t < yi j . But, the matrix X — tS still has zero set Z, so 
our hypothesis on S would be contradicted.) Thus, this part of the integral 
becomes 




(min fe (Y) + ry iojo )e- A < Y+ y™s)d yiojo dY. 



This can be computed as the sum of two integrals in the obvious way. The 
first is di j E(A, Z\J{(io,jo)})/(A'S) and the second is ra,i j /{A-Sf. When 
we sum these expressions over all (zo,jo) £ S we obtain (|T8|). 

When k = m = n < 4, it is easy to see that, when we are not in the 
base case, there always exists a set S of positions in X and disjoint from Z 
such that every minimum ^-assignment of X meets S in the same number of 
positions. 

To illustrate the method we now discuss the case k = m = n = 4. First 
note that if any row or column of X has no fixed zeroes, then we can take 
that row or column to be the set S. So we can suppose that every row or 
column has at least one fixed zero. 

Now suppose that there is a 3 x 3 submatrix S of X that has no fixed 
zeroes. Without loss of generality, we may take this to be the upper left 3x3 
submatrix, so the matrix X has the form 

X = 



where * means that entry is positive and . means that nothing is known 
about that entry. Then any minimum 4-assignment must use either two 
or three entries from S. If it uses three, then it must use 244 and some 
entry x^, i,j < 4. But then we can decrease the value of the assignment 
by replacing and X44 with and x^j, both of which are zero. This 



* 

* 

* 
. 



(19) 
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contradicts the minimality of the 4-assignment we started with. Thus, any 
minimum 4-assignment must use exactly two entries from S. 

If every 3x3 submatrix of X has a fixed zero, and we are not in the 
base case, then the Hall marriage theorem implies that there is a 2 x 3 or 
3x2 submatrix S that has no fixed zeroes. Suppose the former, which we 
can take to be the upper left 2x3 submatrix of X. Each of the first three 
columns has at least one fixed zero. The fixed zeros in those columns must 
be in more than one row, since every 3x3 submatrix has a fixed zero. Thus, 
we may assume the matrix X has the form 



X 



* * * 

* * * 
0.. 
. . . 



(20) 



Any minimum 4-assignment must use one or two entries from S. Suppose 
a minimum 4-assignment uses two entries from S. It cannot use xn, since 
then there would be a smaller 4-assignment consisting of xn,X24,x 3 2,and x i3 . 
Similarly it cannot use x^, x%x, or x 2 2- But, it can only use one of X\3 and 
x 2 3, so it must use exactly one entry from S. 

Thus, we can always find a suitable set S to continue the recursive cal- 
culation. 

We have used this method to compute the expected minimum fc-assignment 
for various small cases. For the case k < m = n < 3, this is easily carried 
out by Mathematica and confirms our conjecture. 

When k = m = n = 4 and k = m = n = 5we were not patient enough 
to wait for Mathematica to simplify the complete rational expression, even 
when the rate matrix has rank 1. However, we were able to check that we 
obtained the correct answer for many random choices or r^'s and Cj's. For this 
purpose, we used an ordinary C program, but, instead of using exact rational 
arithmetic, we carried out our calculations modulo a large prime. Even so, 
the evidence seems to be overwhelming that our conjecture is correct in these 
cases. 

It is possible, although somewhat more complicated, to compute the ex- 
pected contribution of a (k — 1) X (k — 1) submatrix to the expected minimum 
assignment of a k x k matrix when 2 < k < 4. In the cases k = 2 and k = 3 
we were able to check directly with Mathematica that Conjecture |3| was valid, 
which proves Conjecture [l] and Conjecture || whenever k < 3. 
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When k = 4 we obtained computational evidence for the validity of Con- 
jecture checking its validity in a large number of random cases modulo a 
prime. This provides confirmation of the other conjectures when k = 4 and 
m and n are arbitrary. 

5 Additional evidence for the main conjec- 
ture 

Let A = ijiCj) as usual and denote E k (A) by E(k, ri, . . . , r m , c 1; . . . , c n ), 
or simply E(k,r,c). Recall from Section |3.1| the notation F(k,r,c) for the 
formula in Conjecture |2|. In this section we will show that E and F share 
several properties. 

Let us consider how E behaves if we let a collection of the ?Ys approach 

0. For simplicity we assume that n, r2, . . . , ri approach 0. Then the random 
matrices will have very large entries in the first I rows. When k < m — 

1, there are assignments which avoid the first I rows, so in the limit that 
ri, r%, ■ ■ ■ , Ti — > 0, the minimum assignment will avoid those rows and become 
equal to E(k, r i+1 , . . . , r m , c). 

Now suppose that k > m — I. Then a /c-assignment must use at least 
k—(m—l) = k + l — moi the first I rows. But, in our limiting case, these rows 
will be very large, so the minimum fc-assignment will use as few as possible, 
or exactly k + 1 — m of them. The contribution of the entries from these rows 
to the minimum fc-assignment will dominate the minimum k- assignment, so 
in the minimum fc-assignment this contribution will be as small as possible. 
In particular, in the limit as r 1; r 2 , . . . , r; — > 0, this part of the minimum k- 
assignment will be E(k + l — m, ri, . . . , r/, c±, . . . , c n ). By Theorem [| we know 
that a set K of k + / — m columns will be used by the part of the assignment 
in the first / rows with probability Pr(c, K). When this happens the expected 
contribution from the remaining rows is E(m — I, rj+i, . . . , r m , c'(K)) where 
by c'(K) denotes the Cj's corresponding to columns not in K. Thus, we 
should have the following 

Theorem 5 When k < m — I, 

lim E(k,r, c) = E(k,r i+1 , . . . ,r m ,c). (21) 
n,...,^— >o 
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When k > m - I > 0, 



lim (E(k, r, c) — E{k + / — m, r 1; . . . , r h c)) 

ri,...,r t ^O 

= Pr (c, K)E{m - I, r l+1 , ...,r m , c'{K)) (22) 

K 

where the sum is over K C [n] such that \K\ = k + I — m. 

Proof: Let Z be a random exponential mxn matrix with all entries of mean 
1. Then define X = Z/A to be the term by term quotient of the random 
matrix Z by the fixed rate matrix A, where a^- = TjCj. Then X is a random 
exponential matrix with rate matrix A. In particular, 

E(k,r,c) = E k (A) = E(mm k (Z/A)). 

Let A u denote the / x n matrix that is comprised of the first / rows of A, 
and let A d denote (m — I) x n matrix consisting of the last m — I rows of A. 
Furthermore, let Z u and Z^ denote random exponential matrices of the same 
corresponding shapes, again with rate 1. 

If k < m — I, then it is easy to see that 

lim mm k (Z/A) = min k (Z d /A d ) 

ri,...,r;^0 

pointwise almost everywhere (i.e., almost surely, as random variables). Since 
min k (Z/A) < mm k (Z d /A d ), and the expectation of mm k (Z d /A d ) is finite, we 
can apply the dominated convergence theorem to show that the limit of the 
expectation is equal to the expectation of the limit. Thus, 

lim E(mm k (Z/A)) = E(mm k (Z d /A d )) 

ri,...,r ; ^0 

and therefore 

lim E(k,r,c) = E(k,r l+1 , . . . ,r m ,c). 

n,...,ri-»0 

Now suppose k > m — I > 0. Consider the nonnegative random variable 
R = min k (Z / A) —min k+ i_ m (Z u / A u ) . (It is nonnegative because the minimum 
fc-assignment of Z/A must use at least k + I — m elements from Z u /A u .) 
Furthermore, 

Ymn k (Z/A) < mm k+ i_ m (Z u /A u ) + m&x m „i(Z d / A d ) 
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where max^X denotes the maximum ^-assignment of a matrix X. Then 
the second summand is a nonnegative random variable of finite expectation, 
independent of ri, . . . , ri, dominating R. Let xk denote the random variable 
that is 1 or depending upon whether the minimum {k + I — m)-assignment 
of the matrix Z u /A u uses precisely the columns from the set K or does not. 
Then, by Theorem |], we know that E(\k) — Pr(c, K), independent of r. 
Let CK{Z d /A d ) denote the submatrix of Z d /A d obtained when the columns 
indexed by K are removed. We can then show 

lim R - \^ XKmm m _i(C K (Z d / A d )) = 

with convergence pointwise almost everywhere, where the sum is over K C 
[n] such that \K\ = k + / — m. Furthermore, we can bound the finite 
sum by a random variable independent of r 1; . . . ,r^. (The random variable 
™in- m _i(CK(Z d / A d )) will do.) Thus, by the dominated convergence the- 
orem, we can take the limit of the expectations and obtain 



lim E(R) = lim Y^E{ X K^ m -i{C K {Z d / A cl ))) 

= lim Y^E{ XK )E{mxn m ^{C K {Z d /A d ))) 

r\,....ri— >0 ' 
K 

= Pr ( c > K ) E ( m - ^ n +1 , . . . , r m , c'{K)) 

K 

But, E(R) = E(k, r, c) — E{k + I — m, r 1 , . . . , r ; , c), so we are done.D 

We stated Theorem ^ in terms of limits as the first / r's approach 0, in 
order to simplify the notation. However, since E(k,r,c) is symmetric in the 
r's and c's, the analogous results hold for any set of / r's approaching 0. 

From Example p], we have that 

E(l,r h ... ,r m ,d,... ,c n ) = -= — r. (23) 

\2^ii r i)\l^j C j) 

Also when 1 = 1, taking symmetry into account, Theorem |5] reduces to the 
following. 

Suppose that m > 1 and 1 < i < m. Then, if k < m, 

lim E(k, r, c) = E(k, r u . . . , f h . . . , r m , c) (24) 
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while if k = m, 

lim I E{k, r, c) 



n-*0 



Ej c i^( fc - . . ■ , . . . ,r fc ,c x , . . . , Cj, . . . , c n ) 

E; c ; 



(25) 



Proposition 6 There is at most one set of functions 

G(k, . . . , r m , Ci, . . . , c n ), 

each a linear combination of terms of the form ffl), that satisfy the equations 
(Hp, (§3p and /\2^) with E replaced by G. 



Proof: Let H(k,m,n) = H(k,r\,... ,r m ,c%,... ,c n ) denote the difference 
between two sets of functions satisfying the conditions of the proposition. 
We will show by induction on k and m that H(k, m, n) = 0. 

It is clear that H(l,m,n) = for all m and n. Given values of k and 
m, we may suppose that H(k',m',n) = and H(k,m' ,n) = for k' < k, 
m! < m, and arbitrary n. Then, by equations ( |2~3| ) and fl2"5] ) and induction, 
we have, for any i, lim r .^ H(k, m, n) = 0. Suppose that 

H(k,m,n) = J2 ^ ~ . 

where / runs over nonempty subsets of [m] and the /i/s are rational functions 
of the c's. Suppose that hj ^ for some /. Let I be a minimal I such that 
hj ytz o and let t G Jo- Since lim rt ^o H(m, n, k) exists, we must have h\t} = 0, 
so I ^ {t}. Also 

lim H(k, m, n) = ^ — 7 — 

j Et 6 /-{t} r * 

where the sum is over all / strictly containing \t}. Since the terms — are 
linearly independent, we must have hj + hj_{ t y = for all / strictly containing 
{t}. In the case / = Io, this contradicts the minimality of Iq. Thus hj = 
for all / and H(k, m, n) = 0. □ 
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Now we show that the rational functions F(k, r 1? . . . , r m , C\, . . . , c n ) sat- 
isfy the same limit conditions that are proved about E(k, r 1; . . . , r m , ci, . . . , c n ) 
in Theorem [5]. 

In particular the functions F(k, r, c) satisfy the conditions of Proposi- 
tion p. Thus they are the only possible linear combinations of terms of the 
form (H) that could equal E(k,r,c). 

The fact that the F's satisfy all the limit conditions proved about E 
provides additional evidence for Conjecture [2]. 

Theorem 6 When k < m — I, 

lim F(k,r, c) = F(k,r i+1 , . . . ,r m ,c). (26) 

ri,...,n— >0 

When k > m — I > 0, 

lim {F(k, r, c) — F(k + I — m, ri, . . . , 77, c)) 

r 1 ,...,ri^0 

= Pr (c, K)F{m - I, r l+1 , ...,r m , c'{K)) (27) 

K 

where the sum is over K C [n] such that \K\ = k + I — m. 
Proof: We use the alternate form of F(k,r,c) given by 

„ , , . fk — 1— m — n\ 1 

F(*,r,c) = 5:( ifc _ 1 _ J 



, , ^-i-|/|-|/|7(E^n)(E rtJ %)' 

where the binomial coefficient enforces the condition |/| + | J\ < k. 

Now, on the left side of (|2"7]), before passing to the limit, the first term is 

k — 1 — m — n \ 1 



k- 1- III 



where the sum is over subsets / C [m] and J C [n] and the second term is 
y-v / k -i I — hi — ]—/ — ;/ \ I 

i, J 



k + l-m-l-\I\-\J\J (E^/,;<z r *)(E^j9 



where the sum is over subsets / C [/] and JC [n]. By substituting / U {/ + 
1, . . . , m} for / in the second term we get 

y-v / k — ] — ;// — n \ ] 



k-i-\i\-\J\J(E^^(E^jc,) 



33 



where the sum is over subsets I C [m] and J C [n] such that {I + 1, . . . , m} C 
/. Thus, before taking the limit, the left side of (p7|) equals 



A; — 1 — m — n 



k-i-\i\-\J\J (EitinXEjtjCj) 



(28) 



where the sum is over all subsets J C [n] and those subsets / C [m] that do 
not contain {I + 1, . . . , m}. 

The last condition on / implies that n is nonzero if we set r 1; . . . , r; = 
0. Thus, we can obtain the limit on the left of (|2~7D simply by replacing 
ri, . . . ,ri by zero in (|28|). When we do this, the effect is that we combine 
terms with / having a fixed intersection with {I + 1, . . . , m}. Note that this 
intersection is never all of {/ + 1, . . . , m}. 

Suppose then that / is a set strictly contained in {/ + 1, . . . , m}. For each 
i there are (') ways of extending / to a (|/| +i)-element subset of [m] whose 
intersection with {/ + 1, . . . , m} is /. Thus, after taking the limit on the left 
of (P?D, we obtain 

k — 1 — m — n \\ 1 



I, J \i=0 




i) \k - 1 - |/| - \J\ - i) J (Eiti,i>in)(ErfjC 



k + I — 1 — m — n\ 1 

k-i-\i\-\J\) rOtE^c,-) 



where the sum is over proper subsets / C {/ + 1, . . . ,m} and all subsets 
JC [n]. Note that in the case that k < m — I, this expression is precisely 
F(k, ri+i, . . . , r m , c) so we have proved (p6|) . 

We continue with the proof of (|2~TD. We obtain a slightly more convenient 
expression if we replace / by / U [I] in the preceding expression. Then the 
left side becomes 

Ef k + I — 1 — m — n \ 1 
I h 4- j - 1 - i ri - L7i ) rv . ..t-.W . IT^i ^ 29 ^ 



+ /-l-|/|-|J|;(E^r,)(E^ J c J ) 



where the sum is over sets / strictly contained in [m] and containing [I] and 
J C [n). 

Now we turn to the right side of (f27|). 
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The expression (|8|) gives F(m—l, r ;+1 , . . . , r m , c'(K)) as a sum over certain 
subsets I of {I + 1, . . . , m} and J of c'(-fT). But this expression is simpler if 
we replace / by / U [I] and J by J U K. Then the right side can be written 



T3 / „s( k + l-l-m-n\ 



/C[m],_ftTCJC[ri 



where the sum is over proper subsets / C [m] containing [/], subsets C 
J C [n] such that = k + I — m. 

Now we have shown that both the left and right sides of (|27|) are linear 
combinations of the same reciprocal sums l/(Eig/ r «)> so to P rove (0) ^ wm 
suffice to prove that the coefficients of the same reciprocal sums are equal on 
both sides. 

For a given I C [m], the coefficient on the left and right depend only on 
the cardinality of /. We introduce the abbreviations H = k + 1 — 1 — |/| and 
L = k + I — m. Then < L < H < k. After using these abbreviations and 
equating coefficients we are reduced to proving 



JC[n] V \ J \ J ^HJ C 3 KCJC[n] V ^ ' ' ' ^'jt J 3 

where in the sum on the right the subset K must have cardinality L. 
Now use the expression (O) for Pr(c, X) to rewrite the right side of 

as 



(30) 




L-\A\ 




L-n-l\ 1 



H~\J\ JE^jC 



where we still require that \K\ = L. We sum this first over K. In the term 
J2j^ K Cj, the number of K's for which a given Cj occurs depends only on 
whether j belongs to J. Thus, we can rewrite the right side of (BO) as 



E<-i: 



^L-\A\ 

ACJ 



(\J\-\A\\y* , (\J\-\A\-1\ 

I, L-|A| J Z^j^J h ^ \ L-\A\ ) l^j 

E r /,l O ) (E r /./0 




35 



The numerator of the fraction can be rewritten as 




Now, comparing with the left side of (|30|), it suffices to show that the 
inner sum m (§!]) equals (-1) L "I A I ( ) when \A\ < L and otherwise, 
and the inner sum in (^) equals 1. 

It is easy to see that the inner sum of fl3lD equals when \A\ = L, since 
the first binomial coefficient has a negative lower term in that case. When 
\A\ < L, because \J\ > L, the inner sum is over J strictly larger than A. 
Collecting terms according to the cardinality j of J, we obtain 
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H 



E 

j=\A\+l 
H-\A 



n - \A\\ (j - \A\ - l\ (L - n - 1 
j-\A\)\L-\A\-l)[ H-j 



E 



i=l 



n - \A\ 



i - 1 



L-n-1 



L — \ A\ — 1 J \H — \A\ —i 
n- \A\\ ( 0-1 \f L-n-1 

o ){l-\a\-i){h-\a\-o 

n- \A\\ ( i-1 \ f L-n-1 
i ) \L-\A\ -l) \H-\A\ -i 

-1 \fL-n-l\ fn-L + H -\A\ 
L- \A\ - l) V H- \A\ ) + V H-\A\ 
L-n-1 



H-\A\ 



i=0 



\L-\A\ 



H-\A\ 



L-l-H 



where the third equality holds by substituting into the identity (0, p. 16) 
/ m\ / n\ (n + i\ f m — n + q\ fn — m + p\ 

\pJ\qJ ^\p + qJ\ i )\ P-* / 

and the fourth equality holds because L — 1 — H < 0. 

In the inner sum of ([32|) we can collect terms according to the cardinality 
a of A. Recalling that A must be contained in J, we obtain 

This proves Theorem ^. □ 

Finally we prove that E and F have another property in common. 

Theorem 7 i?oin r, c) and r, c) are monotonically decreasing func- 
tions of r and c. In particular, if r x < r[, then 

E(k, ri, r 2) . . . , r m , c) > E(k, r[, r 2) . . . , r m , c) 
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and 



F(k, ri, r 2 , . . . , r m , c) > r' 1; r 2 , . . . , r m , c). 
Furthermore, F is differentiable to any degree £ > in each ri and Cj, and 

(-i)'ff > 0- 

Proof: Without loss of generality, because of symmetry, we can restrict our- 
selves to considering the behavior of E and F as functions of r-y. 

Recall from the proof of Theorem [| that E(k, r, c) is the expectation of 
the random variable mm k (Z/A), where Z is an m x n-matrix- valued random 
variable with exponentially distributed independent entries of mean 1, where 
A = (rjCj) is the rank 1 rate matrix, and where Z/A denotes the element 
by element quotient. Suppose r\ < r[ and = r[ for % = 2, . . . ,m. Define 
A' = (r'iCj). Then A < A' term by term, so Z/A > Z/A', and mmk(Z/A) > 
mxn k (Z/A'). Hence, 

E(k,r,c) = E(mm k (Z/A)) > E{mm k (Z/A')) = E(k,r',c). 

Since there is a nonzero probability that the minimum /^-assignment of Z/A 
uses the first row, the inequality is actually strict. 

We start with the following formula for F, from (|9]) and flTU|): 



|/'|+|j'|<fc \icr Z^g/ i ) \jcji 



For £ > 1, J' C [m], and J' C [n], we define the functions 
/(£,r,/') = E(- 1 ) |/ ' h|/| 7v^I 



ICI> 



g{ 



<c,j>) = 5>i)i™_4- 



so that 



F(fc,r,c)= ]T /(l,r,/'Ml,c,J'). (33) 

|/'|+|j'|<As 
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First we prove that f(£, r, I') > 0. We define the partial sum R = Ylitp r «- 
Then 

t i-i exp(f(-J2- J2 rt))dt 

icp iei 

t t-i exp(t(-5»)cft 

ICP i0 

POO 

= E(- 1 )' / ' h|/| / t^eMt(-J2r i ))dt 

ICP i0 

= ( ,-i) l £(-i ) M-HH_!_ 

= (£-l)\f(£,r,I'), 

so f(£,r,I') > as claimed. The proof that g(£,c, J') > is essentially the 
same. 

To finish the proof of the theorem it suffices to show that for i > 1, 
(-A)V(*,r,c)= ^ + WMW), (34) 

1 |J'| + | J'\=k-1, 10' 



because all terms on the right hand side are positive. First, rewrite (|33| ) as 
follows 

F(k,r,c) = (f(l,r,I') + f(l,r,I f U{l}))g(l,c,J') 

\P\+\J'\<k-l,l<£P 

+ f{±,r,I')g{l,c,J'). 

\I'\ + \J'\=k-l,10' 

The functions f(l,r,F) + f(l, r, I' U {1}) and g(l, c, J') are independent of 
ri, so those partial derivatives vanish. Meanwhile, if % then all the 

denominators in /(l,r, i 7 ) involve r\, so we have 

(-J^yf(l,r,I')=£\f(£+l,r,I'). 
ori 
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This proves (|3~4] ) and Theorem [7|. □ 

The functions E(k, r, c) and F(k, r, c) share many properties. Both func- 
tions are rational, homogeneous of degree —1, and symmetric in the r^'s and 
the c/s. Both rational functions have denominators that factor into linear 
factors which are either sums of r^'s or sums of c/s. Both are positive, and 
both are monotonically decreasing in each and Cj. Finally, E and F share 
various limit properties with or Cj tending to 0. We could also consider 
limits as m or n go to oo. However, it is conceivable that the properties we 
have already found are sufficient to guarantee that such a function is unique. 
In any case, one plan to prove E = F would be to extend the list of common 
properties until equality is forced. 
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